Higher global gross primary productivity under future climate with more advanced representations of photosynthesis.

Gross primary productivity (GPP) is the key determinant of land carbon uptake, but its representation in terrestrial biosphere models (TBMs) does not reflect our latest physiological understanding. We implemented three empirically well supported but often omitted mechanisms into the TBM CABLE-POP: photosynthetic temperature acclimation, explicit mesophyll conductance, and photosynthetic optimization through redistribution of leaf nitrogen. We used the RCP8.5 climate scenario to conduct factorial model simulations characterizing the individual and combined effects of the three mechanisms on projections of GPP. Simulated global GPP increased more strongly (up to 20% by 2070-2099) in more comprehensive representations of photosynthesis compared to the model lacking the three mechanisms. The experiments revealed non-additive interactions among the mechanisms as combined effects were stronger than the sum of the individual effects. The modeled responses are explained by changes in the photosynthetic sensitivity to temperature and CO2 caused by the added mechanisms. Our results suggest that current TBMs underestimate GPP responses to future CO2 and climate conditions.

Detection and attribution of an anomaly in terrestrial photosynthesis in Europe during the COVID-19 lockdown.

Carbon dioxide (CO2) uptake by plant photosynthesis, referred to as gross primary production (GPP) at the ecosystem level, is sensitive to environmental factors, including pollutant exposure, pollutant uptake, and changes in the scattering of solar shortwave irradiance (SWin) - the energy source for photosynthesis. The 2020 spring lockdown due to COVID-19 resulted in improved air quality and atmospheric transparency, providing a unique opportunity to assess the impact of air pollutants on terrestrial ecosystem functioning. However, detecting these effects can be challenging as GPP is influenced by other meteorological drivers and management practices. Based on data collected from 44 European ecosystem-scale CO2 flux monitoring stations, we observed significant changes in spring GPP at 34 sites during 2020 compared to 2015-2019. Among these, 14 sites showed an increase in GPP associated with higher SWin, 10 sites had lower GPP linked to atmospheric and soil dryness, and seven sites were subjected to management practices. The remaining three sites exhibited varying dynamics, with one experiencing colder and rainier weather resulting in lower GPP, and two showing higher GPP associated with earlier spring melts. Analysis using the regional atmospheric chemical transport model (LOTOS-EUROS) indicated that the ozone (O3) concentration remained relatively unchanged at the research sites, making it unlikely that O3 exposure was the dominant factor driving the primary production anomaly. In contrast, SWin increased by 9.4 % at 36 sites, suggesting enhanced GPP possibly due to reduced aerosol optical depth and cloudiness. Our findings indicate that air pollution and cloudiness may weaken the terrestrial carbon sink by up to 16 %. Accurate and continuous ground-based observations are crucial for detecting and attributing subtle changes in terrestrial ecosystem functioning in response to environmental and anthropogenic drivers.

Contrasting anatomical and biochemical controls on mesophyll conductance across plant functional types.

Mesophyll conductance (gm ) limits photosynthesis by restricting CO2 diffusion between the substomatal cavities and chloroplasts. Although it is known that gm is determined by both leaf anatomical and biochemical traits, their relative contribution across plant functional types (PFTs) is still unclear. We compiled a dataset of gm measurements and concomitant leaf traits in unstressed plants comprising 563 studies and 617 species from all major PFTs. We investigated to what extent gm limits photosynthesis across PFTs, how gm relates to structural, anatomical, biochemical, and physiological leaf properties, and whether these relationships differ among PFTs. We found that gm imposes a significant limitation to photosynthesis in all C3 PFTs, ranging from 10-30% in most herbaceous annuals to 25-50% in woody evergreens. Anatomical leaf traits explained a significant proportion of the variation in gm (R2 > 0.3) in all PFTs except annual herbs, in which gm is more strongly related to biochemical factors associated with leaf nitrogen and potassium content. Our results underline the need to elucidate mechanisms underlying the global variability of gm . We emphasise the underestimated potential of gm for improving photosynthesis in crops and identify modifications in leaf biochemistry as the most promising pathway for increasing gm in these species.

Do 2 H and 18 O in leaf water reflect environmental drivers differently?

We compiled hydrogen and oxygen stable isotope compositions (δ2 H and δ18 O) of leaf water from multiple biomes to examine variations with environmental drivers. Leaf water δ2 H was more closely correlated with δ2 H of xylem water or atmospheric vapour, whereas leaf water δ18 O was more closely correlated with air relative humidity. This resulted from the larger proportional range for δ2 H of meteoric waters relative to the extent of leaf water evaporative enrichment compared with δ18 O. We next expressed leaf water as isotopic enrichment above xylem water (Δ2 H and Δ18 O) to remove the impact of xylem water isotopic variation. For Δ2 H, leaf water still correlated with atmospheric vapour, whereas Δ18 O showed no such correlation. This was explained by covariance between air relative humidity and the Δ18 O of atmospheric vapour. This is consistent with a previously observed diurnal correlation between air relative humidity and the deuterium excess of atmospheric vapour across a range of ecosystems. We conclude that 2 H and 18 O in leaf water do indeed reflect the balance of environmental drivers differently; our results have implications for understanding isotopic effects associated with water cycling in terrestrial ecosystems and for inferring environmental change from isotopic biomarkers that act as proxies for leaf water.

Uncovering the critical soil moisture thresholds of plant water stress for European ecosystems.

Understanding the critical soil moisture (SM) threshold (θcrit ) of plant water stress and land surface energy partitioning is a basis to evaluate drought impacts and improve models for predicting future ecosystem condition and climate. Quantifying the θcrit across biomes and climates is challenging because observations of surface energy fluxes and SM remain sparse. Here, we used the latest database of eddy covariance measurements to estimate θcrit across Europe by evaluating evaporative fraction (EF)-SM relationships and investigating the covariance between vapor pressure deficit (VPD) and gross primary production (GPP) during SM dry-down periods. We found that the θcrit and soil matric potential threshold in Europe are 16.5% and -0.7 MPa, respectively. Surface energy partitioning characteristics varied among different vegetation types; EF in savannas had the highest sensitivities to SM in water-limited stage, and the lowest in forests. The sign of the covariance between daily VPD and GPP consistently changed from positive to negative during dry-down across all sites when EF shifted from relatively high to low values. This sign of the covariance changed after longer period of SM decline in forests than in grasslands and savannas. Estimated θcrit from the VPD-GPP covariance method match well with the EF-SM method, showing this covariance method can be used to detect the θcrit . We further found that soil texture dominates the spatial variability of θcrit while shortwave radiation and VPD are the major drivers in determining the spatial pattern of EF sensitivities. Our results highlight for the first time that the sign change of the covariance between daily VPD and GPP can be used as an indicator of how ecosystems transition from energy to SM limitation. We also characterized the corresponding θcrit and its drivers across diverse ecosystems in Europe, an essential variable to improve the representation of water stress in land surface models.

Shifts in the Abundances of Saprotrophic and Ectomycorrhizal Fungi With Altered Leaf Litter Inputs.

Ectomycorrhizal (EcM) and saprotrophic fungi interact in the breakdown of organic matter, but the mechanisms underlying the EcM role on organic matter decomposition are not totally clear. We hypothesized that the ecological relations between EcM and saprotroph fungi are modulated by resources availability and accessibility, determining decomposition rates. We manipulated the amount of leaf litter inputs (No-Litter, Control Litter, Doubled Litter) on Trenched (root exclusion) and Non-Trenched plots (with roots) in a temperate deciduous forest of EcM-associated trees. Resultant shifts in soil fungal communities were determined by phospholipid fatty acids and DNA sequencing after 3 years, and CO2 fluxes were measured throughout this period. Different levels of leaf litter inputs generated a gradient of organic substrate availability and accessibility, altering the composition and ecological relations between EcM and saprotroph fungal communities. EcM fungi dominated at low levels of fresh organic substrates and lower organic matter quality, where short-distances exploration types seem to be better competitors, whereas saprotrophs and longer exploration types of EcM fungi tended to dominate at high levels of leaf litter inputs, where labile organic substrates were easily accessible. We were, however, not able to detect unequivocal signs of competition between these fungal groups for common resources. These results point to the relevance of substrate quality and availability as key factors determining the role of EcM and saprotroph fungi on litter and soil organic matter decay and represent a path forward on the capacity of organic matter decomposition of different exploration types of EcM fungi.

Growth resistance and resilience of mixed silver fir and Norway spruce forests in central Europe: Contrasting responses to mild and severe droughts.

Extreme droughts are expected to increase in frequency and severity in many regions of the world, threatening multiple ecosystem services provided by forests. Effective strategies to adapt forests to such droughts require comprehensive information on the effects and importance of the factors influencing forest resistance and resilience. We used a unique combination of inventory and dendrochronological data from a long-term (>30 years) silvicultural experiment in mixed silver fir and Norway spruce mountain forests along a temperature and precipitation gradient in southwestern Germany. We aimed at examining the mechanisms and forest stand characteristics underpinning the resistance and resilience to past mild and severe droughts. We found that (i) fir benefited from mild droughts and showed higher resistance (i.e., lower growth loss during drought) and resilience (i.e., faster return to pre-drought growth levels) than spruce to all droughts; (ii) species identity determined mild drought responses while species interactions and management-related factors strongly influenced the responses to severe droughts; (iii) intraspecific and interspecific interactions had contrasting effects on the two species, with spruce being less resistant to severe droughts when exposed to interaction with fir and beech; (iv) higher values of residual stand basal area following thinning were associated with lower resistance and resilience to severe droughts; and (v) larger trees were resilient to mild drought events but highly vulnerable to severe droughts. Our study provides an analytical approach for examining the effects of different factors on individual tree- and stand-level drought response. The forests investigated here were to a certain extent resilient to mild droughts, and even benefited from such conditions, but were strongly affected by severe droughts. Lastly, negative effects of severe droughts can be reduced through modifying species composition, tree size distribution and stand density in mixed silver fir-Norway spruce forests.

Sensitivity of gross primary productivity to climatic drivers during the summer drought of 2018 in Europe.

In summer 2018, Europe experienced a record drought, but it remains unknown how the drought affected ecosystem carbon dynamics. Using observations from 34 eddy covariance sites in different biomes across Europe, we studied the sensitivity of gross primary productivity (GPP) to environmental drivers during the summer drought of 2018 versus the reference summer of 2016. We found a greater drought-induced decline of summer GPP in grasslands (-38%) than in forests (-10%), which coincided with reduced evapotranspiration and soil water content (SWC). As compared to the 'normal year' of 2016, GPP in different ecosystems exhibited more negative sensitivity to summer air temperature (Ta) but stronger positive sensitivity to SWC during summer drought in 2018, that is, a stronger reduction of GPP with soil moisture deficit. We found larger negative effects of Ta and vapour pressure deficit (VPD) but a lower positive effect of photosynthetic photon flux density on GPP in 2018 compared to 2016, which contributed to reduced summer GPP in 2018. Our results demonstrate that high temperature-induced increases in VPD and decreases in SWC aggravated drought impacts on GPP. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

Altered energy partitioning across terrestrial ecosystems in the European drought year 2018.

Drought and heat events, such as the 2018 European drought, interact with the exchange of energy between the land surface and the atmosphere, potentially affecting albedo, sensible and latent heat fluxes, as well as CO2 exchange. Each of these quantities may aggravate or mitigate the drought, heat, their side effects on productivity, water scarcity and global warming. We used measurements of 56 eddy covariance sites across Europe to examine the response of fluxes to extreme drought prevailing most of the year 2018 and how the response differed across various ecosystem types (forests, grasslands, croplands and peatlands). Each component of the surface radiation and energy balance observed in 2018 was compared to available data per site during a reference period 2004-2017. Based on anomalies in precipitation and reference evapotranspiration, we classified 46 sites as drought affected. These received on average 9% more solar radiation and released 32% more sensible heat to the atmosphere compared to the mean of the reference period. In general, drought decreased net CO2 uptake by 17.8%, but did not significantly change net evapotranspiration. The response of these fluxes differed characteristically between ecosystems; in particular, the general increase in the evaporative index was strongest in peatlands and weakest in croplands. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

Borehole Equilibration: Testing a New Method to Monitor the Isotopic Composition of Tree Xylem Water in situ.

Forest water use has been difficult to quantify. One promising approach is to measure the isotopic composition of plant water, e.g., the transpired water vapor or xylem water. Because different water sources, e.g., groundwater versus shallow soil water, often show different isotopic signatures, isotopes can be used to investigate the depths from which plants take up their water and how this changes over time. Traditionally such measurements have relied on the extraction of wood samples, which provide limited time resolution at great expense, and risk possible artifacts. Utilizing a borehole drilled through a tree's stem, we propose a new method based on the notion that water vapor in a slow-moving airstream approaches isotopic equilibration with the much greater mass of liquid water in the xylem. We present two empirical data sets showing that the method can work in practice. We then present a theoretical model estimating equilibration times and exploring the limits at which the approach will fail. The method provides a simple, cheap, and accurate means of continuously estimating the isotopic composition of the source water for transpiration.

Higher than expected CO2 fertilization inferred from leaf to global observations.

Several lines of evidence point to an increase in the activity of the terrestrial biosphere over recent decades, impacting the global net land carbon sink (NLS) and its control on the growth of atmospheric carbon dioxide (ca ). Global terrestrial gross primary production (GPP)-the rate of carbon fixation by photosynthesis-is estimated to have risen by (31 ± 5)% since 1900, but the relative contributions of different putative drivers to this increase are not well known. Here we identify the rising atmospheric CO2 concentration as the dominant driver. We reconcile leaf-level and global atmospheric constraints on trends in modeled biospheric activity to reveal a global CO2 fertilization effect on photosynthesis of 30% since 1900, or 47% for a doubling of ca above the pre-industrial level. Our historic value is nearly twice as high as current estimates (17 ± 4)% that do not use the full range of available constraints. Consequently, under a future low-emission scenario, we project a land carbon sink (174 PgC, 2006-2099) that is 57 PgC larger than if a lower CO2 fertilization effect comparable with current estimates is assumed. These findings suggest a larger beneficial role of the land carbon sink in modulating future excess anthropogenic CO2 consistent with the target of the Paris Agreement to stay below 2°C warming, and underscore the importance of preserving terrestrial carbon sinks.

Robust Response of Terrestrial Plants to Rising CO2.

Human-caused CO2 emissions over the past century have caused the climate of the Earth to warm and have directly impacted on the functioning of terrestrial plants. We examine the global response of terrestrial gross primary production (GPP) to the historic change in atmospheric CO2. The GPP of the terrestrial biosphere has increased steadily, keeping pace remarkably in proportion to the rise in atmospheric CO2. Water-use efficiency, namely the ratio of CO2 uptake by photosynthesis to water loss by transpiration, has increased as a direct leaf-level effect of rising CO2. This has allowed an increase in global leaf area, which has conspired with stimulation of photosynthesis per unit leaf area to produce a maximal response of the terrestrial biosphere to rising atmospheric CO2 and contemporary climate change.

Air temperature optima of vegetation productivity across global biomes.

The global distribution of the optimum air temperature for ecosystem-level gross primary productivity ([Formula: see text]) is poorly understood, despite its importance for ecosystem carbon uptake under future warming. We provide empirical evidence for the existence of such an optimum, using measurements of in situ eddy covariance and satellite-derived proxies, and report its global distribution. [Formula: see text] is consistently lower than the physiological optimum temperature of leaf-level photosynthetic capacity, which typically exceeds 30 °C. The global average [Formula: see text] is estimated to be 23 ± 6 °C, with warmer regions having higher [Formula: see text] values than colder regions. In tropical forests in particular, [Formula: see text] is close to growing-season air temperature and is projected to fall below it under all scenarios of future climate, suggesting a limited safe operating space for these ecosystems under future warming.

Quantifying differences in water and carbon cycling between paddy and rainfed rice (Oryza sativa L.) by flux partitioning.

Agricultural crops play an important role in the global carbon and water cycle. Global climate change scenarios predict enhanced water scarcity and altered precipitation pattern in many parts of the world. Hence, a mechanistic understanding of water fluxes, productivity and water use efficiency of cultivated crops is of major importance, i.e. to adapt management practices. We compared water and carbon fluxes of paddy and rainfed rice by canopy scale gas exchange measurements, crop growth, daily evapotranspiration, transpiration and carbon flux modeling. Throughout a monsoon rice growing season, soil evaporation in paddy rice contributed strongly to evapotranspiration (96.6% to 43.3% from initial growth to fully developed canopy and amounted to 57.9% of total water losses over the growing seasons. Evaporation of rainfed rice was significantly lower (by 65% on average) particularly before canopy closure. Water use efficiency (WUE) was significantly higher in rainfed rice both from an agronomic (WUEagro, i.e. grain yield per evapotranspiration) and ecosystem (WUEeco, i.e. gross primary production per evapotranspiration) perspective. However, our results also show that higher WUE in rainfed rice comes at the expense of higher respiration losses compared to paddy rice (26% higher on average). Hence, suggestions on water management depend on the regional water availability (i.e. Mediterranean vs. Monsoon climate) and the balance between higher respiratory losses versus a potential reduction in CH4 and other greenhouse gas emissions. Our results suggest that a shift from rainfed/unsaturated soil to waterlogged paddy conditions after closure of the rice canopy might be a good compromise towards a sustainable use of water while preserving grain yield, particularly for water-limited production areas.

Validating MODIS and Sentinel-2 NDVI Products at a Temperate Deciduous Forest Site Using Two Independent Ground-Based Sensors.

Quantifying the accuracy of remote sensing products is a timely endeavor given the rapid increase in Earth observation missions. A validation site for Sentinel-2 products was hence established in central Germany. Automatic multispectral and hyperspectral sensor systems were installed in parallel with an existing eddy covariance flux tower, providing spectral information of the vegetation present at high temporal resolution. Normalized Difference Vegetation Index (NDVI) values from ground-based hyperspectral and multispectral sensors were compared with NDVI products derived from Sentinel-2A and Moderate-resolution Imaging Spectroradiometer (MODIS). The influence of different spatial and temporal resolutions was assessed. High correlations and similar phenological patterns between in situ and satellite-based NDVI time series demonstrated the reliability of satellite-based phenological metrics. Sentinel-2-derived metrics showed better agreement with in situ measurements than MODIS-derived metrics. Dynamic filtering with the best index slope extraction algorithm was nevertheless beneficial for Sentinel-2 NDVI time series despite the availability of quality information from the atmospheric correction procedure.

Competing influences of anthropogenic warming, ENSO, and plant physiology on future terrestrial aridity.

The 2011-2016 Californian drought illustrates that drought-prone areas do not always experience relief once a favorable phase of El Niño-Southern Oscillation (ENSO) returns. In the 21st century, such an expectation is unrealistic in regions where global warming induces an increase in terrestrial aridity larger than the aridity changes driven by ENSO variability. This premise is also flawed in areas where precipitation supply cannot offset the global warming-induced increased evaporative demand. Here, atmosphere-only experiments are analyzed to identify land regions in which aridity is currently sensitive to ENSO, and where projected future changes in mean aridity exceed the range caused by ENSO variability. Insights into the drivers of these aridity changes are obtained in simulations with incremental addition of three different factors to current climate: ocean warming, vegetation response to elevated CO2 levels, and intensified CO2 radiative forcing. The effect of ocean warming overwhelms the range of ENSO-driven temperature variability worldwide, increasing potential evapotranspiration (PET) in most ENSO-sensitive regions. Additionally, ~39% of the regions currently sensitive to ENSO receive less precipitation in the future, independent of the ENSO phase. Aridity increases consequently in 67-72% of the ENSO-sensitive area. When both radiative and physiological effects are considered, the area affected by aridity rises to 75-79% when using PET-derived measures of aridity, but declines to 41% when total soil moisture aridity indicator is employed. This reduction mainly occurs because plant stomatal resistance increases under enhanced CO2 concentrations, which results in improved plant water use efficiency, and hence reduced evapotranspiration and soil desiccation. Imposing CO2-invariant stomatal resistance may overestimate future drying in PET-derived indices.

An empirical study of the wound effect on sap flux density measured with thermal dissipation probes.

The insertion of thermal dissipation (TD) sensors on tree stems for sap flux density (SFD) measurements can lead to SFD underestimations due to a wound formation close to the drill hole. However, the wound effect has not been assessed experimentally for this method yet. Here, we propose an empirical approach to investigate the effect of the wound healing on measured sap flux with TD probes. The approach was performed for both, diffuse-porous (Fagus sylvatica (Linnaeus)) and ring-porous (Quercus petraea (Lieblein)) species. Thermal dissipation probes were installed on different dates along the growing season to document the effects of the dynamic wound formation. The trees were cut in autumn and additional sensors were installed in the cut stems, therefore, without potential effects of wound development. A range of water pressures was applied to the stem segments and SFDs were simultaneously measured by TD sensors as well as gravimetrically in the laboratory. The formation of wounds around sensors installed in living tree stems led to underestimation of SFD by 21.4 ± 3 and 47.5 ± 3.8% in beech and oak, respectively. The differences between SFD underestimations of diffuse-porous beech and ring-porous oak were, however, not statistically significant. Sensors with 5-, 11- and 22-week-old wounds also showed no significant differences, which implies that the influence of wound formation on SFD estimates was completed within the first few weeks after perforation. These results were confirmed by time courses of SFD measurements in the field. Field SFD values decreased immediately after sensor installation and reached stable values after ~2 weeks with similar underestimations to the ones observed in the laboratory. We therefore propose a feasible approach to correct directly field observations of SFD for potential underestimations due to the wound effect.

The plumbing of land surface models: is poor performance a result of methodology or data quality?

The PALS Land sUrface Model Benchmarking Evaluation pRoject (PLUMBER) illustrated the value of prescribing a priori performance targets in model intercomparisons. It showed that the performance of turbulent energy flux predictions from different land surface models, at a broad range of flux tower sites using common evaluation metrics, was on average worse than relatively simple empirical models. For sensible heat fluxes, all land surface models were outperformed by a linear regression against downward shortwave radiation. For latent heat flux, all land surface models were outperformed by a regression against downward shortwave, surface air temperature and relative humidity. These results are explored here in greater detail and possible causes are investigated. We examine whether particular metrics or sites unduly influence the collated results, whether results change according to time-scale aggregation and whether a lack of energy conservation in flux tower data gives the empirical models an unfair advantage in the intercomparison. We demonstrate that energy conservation in the observational data is not responsible for these results. We also show that the partitioning between sensible and latent heat fluxes in LSMs, rather than the calculation of available energy, is the cause of the original findings. Finally, we present evidence suggesting that the nature of this partitioning problem is likely shared among all contributing LSMs. While we do not find a single candidate explanation for why land surface models perform poorly relative to empirical benchmarks in PLUMBER, we do exclude multiple possible explanations and provide guidance on where future research should focus.

Sauna, sweat and science - quantifying the proportion of condensation water versus sweat using a stable water isotope ((2)H/(1)H and (18)O/(16)O) tracer experiment.

Most visitors of a sauna appreciate the heat pulse that is perceived when water is poured on the stones of a sauna stove. However, probably only few bathers are aware that this pleasant heat pulse is caused by latent heat being released onto our skin due to condensation of water vapour. In order to quantify the proportion of condensation water versus sweat to dripping water of test persons we conducted sauna experiments using isotopically labelled (δ(18)O and δ(2)H) thrown water as tracer. This allows differentiating between 'pure sweat' and 'condensation water'. Two ways of isotope mass balance calculations were applied and yielded similar results for both water isotopes. Accordingly, condensation contributed considerably to dripping water with mean proportions of 52 ± 12 and 54 ± 7% in a sauna experiment in winter semester 2011/12 and 30 ± 13 and 33 ± 6% in a sauna experiment in winter semester 2012/13, respectively, depending on the way of calculating the isotope mass balance. It can be concluded from the results of our dual isotope labelling sauna experiment that it is not all about sweat in the sauna.

Stable oxygen isotope and flux partitioning demonstrates understory of an oak savanna contributes up to half of ecosystem carbon and water exchange.

Semi-arid ecosystems contribute about 40% to global net primary production (GPP) even though water is a major factor limiting carbon uptake. Evapotranspiration (ET) accounts for up to 95% of the water loss and in addition, vegetation can also mitigate drought effects by altering soil water distribution. Hence, partitioning of carbon and water fluxes between the soil and vegetation components is crucial to gain mechanistic understanding of vegetation effects on carbon and water cycling. However, the possible impact of herbaceous vegetation in savanna type ecosystems is often overlooked. Therefore, we aimed at quantifying understory vegetation effects on the water balance and productivity of a Mediterranean oak savanna. ET and net ecosystem CO2 exchange (NEE) were partitioned based on flux and stable oxygen isotope measurements and also rain infiltration was estimated. The understory vegetation contributed importantly to total ecosystem ET and GPP with a maximum of 43 and 51%, respectively. It reached water-use efficiencies (WUE; ratio of carbon gain by water loss) similar to cork-oak trees. The understory vegetation inhibited soil evaporation (E) and, although E was large during wet periods, it did not diminish WUE during water-limited times. The understory strongly increased soil water infiltration, specifically following major rain events. At the same time, the understory itself was vulnerable to drought, which led to an earlier senescence of the understory growing under trees as compared to open areas, due to competition for water. Thus, beneficial understory effects are dominant and contribute to the resilience of this ecosystem. At the same time the vulnerability of the understory to drought suggests that future climate change scenarios for the Mediterranean basin threaten understory development. This in turn will very likely diminish beneficial understory effects like infiltration and ground water recharge and therefore ecosystem resilience to drought.